Force induced conformational transition in a system of interacting 
stiff polymer: Application to unfolding 
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Abstract 

We consider a stiff polymer chain in poor solvent and apply a force at one end of the chain. We find 
^ | that by varying the stiffness parameter, polymer undergoes a transition from the globule state to the 
folded like state. The conformation of folded state mimics the /3-sheet as seen in titin molecule. Using 

0- 

exact enumeration technique, we study the extension-force and force-temperature diagrams of such a 
system. Force-temperature diagram shows the re-entrance behaviour for flexible chain. However, for 
stiff chain this re-entrance behaviour is absent and there is an enhancement in ^-temperature with the 
rise of stiffness. We further propose that the internal information about the frozen structure of polymer 

>. 

can be read from the distribution of end-to-end distance which shows saw-tooth like behaviour. 
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Conformational and structural properties of bio-polymers play an important role in the biolog- 
ical phenomena. Protein folding is a phenomenon which is associated with the primary structure 
of the chain. Their stability or resistance to unfolding are investigated either by varying pH of 
the solution or by varying the temperature Many proteins unfold due to the force instead 
of change in the chemical environment or temperature. For example, the effect of stress on the 
titin molecule is a force induced transition from the native state to the unfolded state [2j. The 
presence of strong hysteresis together with sudden jump seen in the force vs extension curve 
suggests that the unfolding is a first order transition [3]. Theoretically, this type of transition 
has been studied in the framework of statistical mechanics considering a flexible polymer chain 
in a poor solvent under the influence of external force Q], and for DNA type polymers under 
an unzipping force However, stiffness plays an important role as is seen in worm like chain 
(WLC) model 3, S], which has not been incorporated in the lattice model so far to explain the 
phase diagram of such molecules under the influence of the external force. 

Recently the phase diagram of homo semi-flexible polymer chain with zero force has been 
studied by mean field theory and Monte Carlo simulation Ml which shows three distinct 
phases namely : (i) an open coil phase at high temperature, (ii) a molten globule collapse at low 
temperature and low stiffness, and (iii) a 'frozen' or 'folded' state at low temperature and large 
stiffness. Our aim in this paper is to describe the effect of external force applied at one end of 
the stiff polymer chain in poor solvent. We consider self-attracting self-avoiding walks (SAWs) 
on a square lattice One end of the chain is subjected to an external force while the other 
end is kept fixed. The stretching energy E s arising due to the applied force F is taken as 

E s = -F.x (1) 

where x is the x-component of end-to-end distance (| x\ — xn |)- Stiffness in the chain is in- 
troduced by associating an energy barrier e& with every 'turn' (or bend) of the walk 9fl. We 
associate a negative energy e u for each non-bonded occupied nearest neighbour pairs. The par- 
tition function of such a system may be written as 

Z N = C N (N b ,N u ,\x\)b N "u N ^ (2) 

(N b ,N u ,\x\) 

CN(Nb, N u , I x |) is the total number of SAWs [10] of N steps having Nf, turns (bends) and N u 
nearest neighbour pairs respectively, u is the Boltzmann weight for the force which is defined 
as exp[/3(F.x)], where x is the unit vector along the x-axis. [3 is defined as where k is the 
Boltzmann constant and T is the temperature, b = exp(—/3eb) and u = exp(—/3e u ) are Boltzmann 
weights of bending and nearest neighbour interaction respectively. 
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FIG. 1: Phase diagram shows the variation of b (Boltzmann weight of bending) with u (Boltzmann 
weight of nearest neighbour interaction) of a linear semi-flexible polymer chain in two dimensional 
space. The collapsed state includes both globule and folded state. The open circle represents transition 
from swollen state to the collapsed state and the filled circle represents transition from the folded state 
to the globule state. 

Since the polymer is on the lattice where continuous overall rotations are not possible, we 
assume that under the applied force there is an alignment of protein along the force direction with 
zero torque. This is equivalent to the assumption that the relaxation for the overall rotational 
degrees of freedom is much shorter than that of the structural relaxation that is responsible for 
unfolding or folding process. Apart from logarithmic factors, the low temperature limit of the 
rotationally averaged contribution to the the free energy due to tension obtained in the context 
of rubber elasticity ll| coincides with Eq. (1). 

In all the single molecule experiments, a finite size of the chain is used and the fact that no 
true phase transition can occur in a single macromolecule, we calculate "state diagram" instead 
of "phase diagram" [12l |. We obtained Cjv(iV&, N u , \ x |) up to N < 30 steps walk on a square 
lattice by the exact enumeration method and use it to calculate Z^(b,u,uj). The boundaries 
in the state diagram may be found for F = from the maxima of fluctuation of N u (= d< I^ u> ) 



or fluctuation of iVt,(= g e )■ The sudden change in the non-bonded nearest neighbors N u 
indicates a phase transition; hence maxima in the derivative of N u with respect to e u indicates 
the location of a phase transition. The swollen to collapsed transition line (which is also called 
the (9-line) is obtained from the peak value of d< ^ u> for fixed value of b (shown in Fig. 1 by open 
circle). Folded to globule transition line (shown in Fig. 1 by filled circle) has been obtained from 



the peak value of 



d<N b > 



at fixed u. 



There are three states marked by swollen, frozen (or folded) and globule as shown in Fig. 1. 
For b = 1, we restore the value of u c = 1.93 using extrapolation scheme |l3j for flexible polymer 
chain, that is in good agreement with the value found by Foster [ljj] and the Monte Carlo 
simulation result (1.94 ± 0.005) 151 ] . However, we shall confine ourselves here to the constant 
force ensemble to get the exact boundaries for finite N. 

From Fig. 1, we find that the #-line bends to left as stiffness increases (b decreases) indicating 
the enhancement of ^-temperature. Intuitively one might think that stiffness favors the extended 
state but the fact is that stiffness favors a folded state in a poor solvent at low temperature. It 
is obvious that on a square lattice hairpin like structures (Fig. 2(a)) minimize the number of 
folds and maximize the number of nearest neighbors. Due to this increase in nearest neighbours, 
there is a decrease in u (i.e. rise in ^-temperature). Such trend in #-line has also been observed 

Q 

by extensive Monte Carlo simulation by Bastolla and Grassberger [2J. For large N, resulting 
configurations will give rise frozen or folded like state and the partition function will be then 
dominated by configurations similar to Fig. 2(b). Since the number of such configurations are 
very small, entropy associated with folded state will be quite low. The length of the fold will be 
of the order of N 1 / 2 and depends on (/3e&). 

To study the order of transition from the frozen state to the globule state at zero force, we 
calculated the entropy (S) using the following relation of free energy (G): 



G = -kT\nZ N (T) 

_ fdG 
b ~ ~[df 



(3) 
(4) 



We use two different sets of e& at fixed u = 4.0 and plotted the entropy with temperature 
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FIG. 2: Typical conformations of polymer chain at higher value of u and lower value of b in two 
dimensional space: (a) This represents a situation where hairpin structure will be preferred than the 
other conformations; (b) The resulting conformation at large N; (c) The conformation at T = 0. This 
can be mapped by Hamiltonian walk. 
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FIG. 3: Plot of Entropy (S) vs Temperature (T) for different values of stiffness parameter (q,). 



in Fig. 3. In the first set (Fig. 3(a)), we choose e& = —0.08 eV and —0.16 eV, implying 
b = exp(+/3eb) is greater than one, and this corresponds to a flexible chain. For the second set 
(shown in Fig. 3(b)) is positive (b < 1), and this corresponds to the stiff chain. It can be seen 
that for the stiff chain there is a sudden jump in entropy (shown in Fig. 3(b)) corresponding to 
a transition from the folded state to the globule state while it is absent in the case of a flexible 
polymer chain. 

Although in principle the frozen structures (Fig. 2(b)) can be seen everywhere in the "Folded" 
state, numerically it is easier to observe this for smaller b (< 0.2) and larger u (> 1.93). In this 
regime we have numerically computed the Boltzmann contribution of individual conformation 
of polymer chain in the partition function and verified that the maximal contribution comes 
from the structure similar to Fig. 2(b) (like /5-sheet seen in titin molecules). Hence in this 
region, it will be possible to study the effect of force on the unfolding transition by applying 
force at one end of the chain. To do so we set e u = — 1 and Boltzmann constant k — 1 and 
study force-temperature and extension-force curve. It can be shown that in this scale F = ({^), 
T=l/ln( W ) and e b = -(£*). 

The variation of critical force with temperature for unfolding transition where polymer goes 
from the collapsed state to the unfolded state (or extended state) is shown in Fig. 4. Note that 
the collapsed state consists of both the globule and the folded state. The transition line separates 
the collapsed state from the unfolded state. It is not possible to see transition from the folded 
to the globule state because it is not induced by the force. It is interesting to note that for a 
stiff polymer chain, critical force increases monotonically with temperature and becomes almost 
constant at very low temperature. However, in the case of flexible polymer chain, we observe 
a phenomenon of 're-entrance' i.e. the critical force goes through maximum as temperature is 
lowered. For example, one can see that the polymer chain at fixed force (say F = 0.9), is in 



the extended state at low temperature. With the increase in temperature, the chain is found in 
the collapsed state. With the further rise of temperature, it acquires again the conformations of 
the extended (swollen) state. Similar re-entrance behaviour is also found in the transfer matrix 

n 

calculation of the directed walk models of flexible polymer chain [16j. However, if we introduce 
stiffness in the chain, the re-entrance behaviour is found to be suppressed. This is because in 
the globule state, the entropy associated with the flexible chain is very high, while the stretched 
chain has almost zero entropy (i.e polymer forms a rod like shape). This indicates that polymer 
goes from a high entropy state to the low entropy state under the application of force for flexible 
polymer chain. However, in the case of stiff polymer chain, collapsed state has frozen structure 
with almost zero entropy, therefore, no re-entrance is observed in going from the frozen state 
to the extended state. This can be seen using phenomenological argument near T = where 
the conformation of polymer chain in the collapsed state looks like Fig. 2(c) similar to the 
conformations formed by the Hamiltonian walks. Thus, from the principle of balance of energy, 
the free energy of folded state and the free energy of stretched state due to the force (using Eq. 
1) can be equated as [k^ 

-FN = Ne u + 2v / iV(e 6 - e u ) - NTS C (5) 

The second term of right hand side in Eq. (5) is a surface correction term which also includes 
the bending energy. The third term is due to the entropy associated with the collapsed state (S c 
per monomer). Minimization of energy with respect to N and substituting e u = —1, Eq. (5) 
gives 

F C (T) = 1~^S + TS C (6) 



For large N, the second term goes to zero. For finite N = 30 (as we have taken in our study), 
the correction term is important which gives F c = 0.8174. This is in good agreement with the 
value found from Fig. 4 near T = 0. It is to be noted that for a flexible polymer, entropy 
associated with globule is finite and hence there is a positive slope (^), while for stiff chain 
entropy associated with frozen structure is almost zero and hence there is no slope, therefore, we 
do not find any re-entrance behavior in this case. 

The extension vs force (< x > vs F) curve is shown in Fig. 5. From the figure, it is evident 
that for small forces, a polymer chain is in the compact folded state and slightly oriented along 
the force direction. At larger forces, the polymer chain has the conformations similar to the 
extended (swollen) structure. Completely stretched states can be obtained only by applying 
very high forces. A similar behaviour is also seen in the case of flexible polymer chain under 
the influence of external forces Molecular dynamics simulation of protein under the stress 
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FIG. 4: Force- Temperature phase diagram for fixed b. Here collapsed state consists of both the globule 
and the folded state. 
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FIG. 5: Plot of < x > vs F for fixed u and b, Here u = 10.0 and b = 0.2. 

of an external denaturating force acting on a terminal end or on the entire chain [18] shows the 
intermediate stages during the unfolding and the extension-force curves are similar to Fig. 5 
found by us. 

We also study the probability distribution curves P(x) with x for flexible and stiff chains 
defined by 



P(\x\) 



1 



Zpf(b, it, u) t 



J2C N (N b ,N u ,\x\)b 



(7) 



(1*1) 



In Fig. 6, we have shown P(x) for different values of uj and at fixed u = 3.0 corresponding to 
the collapsed state. 

The probability distribution curves have many interesting features. For flexible chain i.e. 
6 = 1, the maxima of P(x) corresponds to the collapsed state at F = (uj = 1.0). For 
higher force, to = 20 both flexible and semi-flexible polymers are found to be in the "rod-like" 
state. However, at intermediate force, the probability distribution curve has "saw-tooth" type 
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FIG. 6: Plot of P{x) vs x for flexible (b = 1.0) and stiff (6 = 0.2)) polymer chain at different u> 
(Boltzmann weight of the force). At intermediate force the saw-tooth behaviour in the probability 
distribution gives signature of unfolding events in the frozen structure. 

of behaviour for the stiff chain while it is continuous for the flexible chain. The x-component 
of end-to-end distribution function gives information about the internal structure of the folded 
state by applying the suitable force. For small forces, the thermal fluctuations are too weak to 
unfold the polymer chain, and it stays in the folded state most of the time. This fact is more or 
less reflected in the structureless distribution function with a well denned peak at the most likely 
value of the end-to-end distance. In contrast, when we apply a force close to the critical force, 
the small thermal fluctuations are sufficient to open up the loops in the /3-sheet (as shown in Fig. 
2(b)). Here, one sees the statistical unfolding events in the structure of the distribution functions 
in the form of the "saw-tooth" kind of behaviour for the stiff chain while it is continuous for the 
flexible chain. The presence of intermediate stages seen during the unfolding in the molecular 



dynamic simulation 



181 ] and the one found by us in the probability distribution curve, suggests 



the experimentalist to "tease" the folded polymer (proteins) by some external force and then 
measure the distribution function. This may reveal some interesting information about the folded 
state. 

In this paper we have studied the complete state diagram of stiff polymer chain under the 
influence of external force. We showed the existence of the folded-like state in stiff homopolymer. 
We have also found that there is an enhancement of ^-temperature (decrease in u) with the rise 
of stiffness parameter e&. The absence of re-entrance in the stiff chain has been explained by 
using the phenomenological argument. We have also tried to explain the unfolding event as seen 
in titin molecule on the basis of the probability distribution curve. 
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